Manipulacion de datos

library(haven)
## Warning: package 'haven' was built under R version 4.0.3
library(tidyverse)
## Warning: package 'tidyverse' was built under R version 4.0.3
## -- Attaching packages --------------------------------------- tidyverse 1.3.0 --
## v ggplot2 3.3.2     v purrr   0.3.4
## v tibble  3.0.4     v dplyr   1.0.2
## v tidyr   1.1.2     v stringr 1.4.0
## v readr   1.4.0     v forcats 0.5.0
## Warning: package 'ggplot2' was built under R version 4.0.3
## Warning: package 'tibble' was built under R version 4.0.3
## Warning: package 'tidyr' was built under R version 4.0.3
## Warning: package 'readr' was built under R version 4.0.3
## Warning: package 'purrr' was built under R version 4.0.3
## Warning: package 'dplyr' was built under R version 4.0.3
## Warning: package 'stringr' was built under R version 4.0.3
## Warning: package 'forcats' was built under R version 4.0.3
## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()
Caracteristicas_y_composicion_del_hogar <- read_sav("C:/Users/julia/Downloads/Encuesta Calidad de Vida/Caracteristicas y composicion del hogar.sav")

Servicios_del_hogar <- read_sav("C:/Users/julia/Downloads/Encuesta Calidad de Vida/Servicios del hogar.sav")

http://microdatos.dane.gov.co/index.php/catalog/607/datafile/F116

Caracteristicas_y_composicion_del_hogar
#extraccion numero de hijos
nhijos <- Caracteristicas_y_composicion_del_hogar %>% 
   filter(P6051==3, SECUENCIA_P == 1) %>%
   group_by(DIRECTORIO) %>%
   summarise(Hijos=n()) 
## `summarise()` ungrouping output (override with `.groups` argument)
barplot(table(nhijos$Hijos), main="Numero de hijos")

caracteristicas_hogar <- Caracteristicas_y_composicion_del_hogar %>% 
   filter(P6051==1, SECUENCIA_P == 1) %>%
   select(DIRECTORIO ,P6020, P6040, P5502, P6071, P5667)
#secuencia encuesta 1 para la primera encuesta que le hicieron al hogar
servicios <- Servicios_del_hogar %>%
   filter(SECUENCIA_ENCUESTA == 1) %>% 
   select(DIRECTORIO, CANT_PERSONAS_HOGAR, I_HOGAR, P5010) 
#unir caracteristicas y nhijos a datos
datos <- merge(x = caracteristicas_hogar, y = nhijos, by = "DIRECTORIO", all.x=TRUE, all.y=TRUE)
#unir datos con servicios 
datos <- merge(x = datos, y = servicios, by = "DIRECTORIO", all.x = TRUE, all.y=TRUE)

#reemplazar NA con 0 hijos
datos <- datos %>% replace_na(list(Hijos = 0))
datos %>%
   filter(CANT_PERSONAS_HOGAR == 0 & Hijos == 0)
plot(datos$CANT_PERSONAS_HOGAR,datos$Hijos)

write.csv(datos,"datos.csv", row.names = TRUE)